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O ; Abstract 

'NT 

S= I We review our results on a mathematical dynamical theory for observables for open many- 
body quantum nonlinear bosonic systems for a very general class of Hamiltonians. We show that 

Mh, non-quadratic (nonlinear) terms in a Hamiltonian provide a singular "quantum" perturbation for 

fl ' observables in some "mesoscopic" region of parameters. In particular, quantum effects result in 



secular terms in the dynamical evolution, that grow in time. We argue that even for open quan- 
tum nonlinear systems in the deep quasi-classical region, these quantum effects can survive after 
decoherence and relaxation processes take place. We demonstrate that these quantum effects in 
open quantum systems can be observed, for example, in the frequency Fourier spectrum of the 
dynamical observables, or in the corresponding spectral density of noise. Estimates are presented 
for Bose-Einstein condensates, low temperature mechanical resonators, and nonlinear optical sys- 
tems prepared in large amplitude coherent states. In particular, we show that for Bose-Einstein 
condensate systems the characteristic time of deviation of quantum dynamics for observables from 
the corresponding classical dynamics coincides with the characteristic time-scale of the well-known 
quantum nonlinear effect of phase diffusion. 

PACS numbers: 03.67.Lx, 75.10.Jm 



I. INTRODUCTION 

Real physical systems are not isolated, they are coupled to external degrees of freedom. 
The classical and quantum dynamics of these open systems are especially complex for non- 
linear systems that exhibit several phenomena, including deviation of quantum dynamics 
from the corresponding classical one, quantum revivals, decoherence, and relaxation. Re- 
cently substantial effort has been devoted to study the open dynamics of nonlinear quantum 
systems, with the aim of understanding the quantum to classical transition in a controlled 
way pj 

Standard mathematical treatments of open quantum nonlinear systems suffer from prob- 
lems arising from the interplay between the nonlinearity and the openness of the system. 
Usually the dynamics of open quantum systems is studied using different mathematical ap- 
proaches, such as the master equation for the reduced density matrix, which is an average 
of the full density matrix over the environment. [2, y, |j], and quasi-probability distribu- 
tions (e.g. the so-called Q-function 5|], the Wigner function 6(], etc). Although all of these 
approaches allow one, in principle, to calculate the time evolution of the average values of 
the dynamical variables of the system, they have significant drawbacks. In particular, these 
distribution functions may not be positively defined; they may be inconsistent for certain 
density matrices; it may be difficult to extract physical information from these distributions, 
especially in the context of quantum nonlinear open systems; in the "deep" quasi-classical 
region of parameters, e = h/ J <C 1 (where h is Planck constant and J is a characteristic 
action of the corresponding classical system) these quasi-probability distributions exhibit 
fast oscillations due to phases like exp(zS'(t)//i), with |5'(t)| ~ J. Therefore, it is difficult 
to separate the physical effects for dynamical observables (requiring an additional multi- 
dimensional integration of quasi-distribution densities) from the effects of errors related to 
a concrete mathematical approach. 

We are approaching these problems using an alternative strategy that starts from a 
mathematical dynamical theory based on exact, linear partial differential equations (PDEs) 
for the observables of open many-body quantum nonlinear bosonic systems governed by a 
very general class of Hamiltonians (see [7|, la, |9| and references therein). The key advantage of 
this method is that it leads to a well-behaved asymptotic theory for open quantum systems 
in the quasi-classical region of parameters. This approach is a generalization to the open 



case of the asymptotic theory for bosonic and spin closed quantum systems [7|, [lO|, lUl , and 
it can be apphed to general open quantum nonlinear bosonic and spin systems for a large 
range of parameters, including the deep quasi-classical region. 

We concentrate our attention on a discussion of the method which can be used to ob- 
serve quantum effects after decoherence and relaxation, in the deep quasi-classical region 
of parameters. We argue that one can use for these purposes a Fourier spectrum of the 
dynamical observables, since its width contains characteristic information of such quantum 
effects. Our observation is based on our first studies 8|, l9| of this new approach to quan- 
tum nonlinear systems interacting with an environment. As will be discussed below, certain 
quantum effects which are presented in the dynamics of these nonlinear systems are robust 
to the influence of the environment, and survive after decoherence and relaxation processes 
take place. In order to observe these effects experimentally it is necessary to have a quasi- 
classical system in certain region of parameters. We call these systems "mesoscopic" , mainly 
because the parameter e should not be too small. In this sense, many quasi-classical sys- 
tems have the drawback that they are either "too classical" (i.e., they have a large J so 
that the quasi-classical parameter e is extremely small), or they interact too strongly with 
the environment, or their effective temperature is so high that quantum effects that we are 
talking about are washed out. Only recently have adequate open nonlinear quasi-classical 
systems become available, including Bose-Einstein condensates (BEG) with large number of 
atoms and thermally well isolated; high frequency cantilevers with large nonlinearities and 
at sufficiently low temperatures; and nonlinear optical systems in high Q resonators, among 
others. We present estimates on the parameter regions where survival of certain quantum 
effects to environment-induced decoherence can be observed in these systems. 

II. DYNAMICS OF QUANTUM OBSERVABLES FOR CLOSED QUANTUM 
NONLINEAR SYSTEMS 

We first consider closed quantum nonlinear systems. As a simple example we take the 



one-dimensional quantum nonlinear oscillator (QNO) described by the Hamiltonian 
(see also an application of this Hamiltonian for the BEG system in Section V) 



0.^ 



M 



Hs = hjJo)a + ^f\^{a)af' , [a\ a\ = 1, (1) 



where a, a^ are the annihilation and creation operators, u is the frequency of hnear os- 
cillations, and /i is a dimensional parameter of nonlinearity. We assume that initially 
the QNO is prepared in a coherent state \a) {a\a) = a\a)). In the classical limit 
(a ^ a,a^ ^ a*, |ap -^ oo,/l|ap = J, the classical action of the linear oscillator) the 
Hamiltonian ([T]) becomes Hd = ujj + fiJ"^. Below we use the following dimensionless nota- 
tion: r = ut, /i = hfi/u, and /id = fiJ/u. The quantum parameter of nonlinearity /i can be 
presented as the product of two parameters, quantum and classical, fl = e/id- The parameter 
/ici characterizes the nonlinearity in the classical nonlinear oscillator (BEC, cantilever, opti- 
cal field, etc) and can be written as /id = {J /2uj) {dud/ dJ) , where Ud = dHd/dJ = u + 2/iJ 
is the classical frequency of nonlinear oscillations. The limit /id ^ 1 corresponds to weak 
nonlinearity, while /id — 1 corresponds to strong nonlinearity. As was mentioned above, e is 
the quasi-classical parameter. Namely, e ~ 1 corresponds to the pure quantum system, and 
e <^ 1 corresponds to the quasi-classical limit, which is the subject of our interest. 

A. Closed partial differential equation for observables 

A closed linear PDE which describes the time evolution of the expectation value of any 
observable of the system can be easily derived when the system is initially populated in a 
coherent state |a) (see [7j and references therein). Namely, for an arbitrary operator function 
/ = f{a,a^), the time-dependent expectation value (observable) of such a function, 

/(a*,a,r) = (a|e^^*/Ve-^^*/V), (2) 

satisfies a PDE of the form 

df/dr = Kf, (3) 

wherein = Kd + e/id-f^q- Here the operator Kd includes only the first order derivatives 
and describes the corresponding classical limit, while the operator K^ includes higher-order 
derivatives and contains the quantum effects. For the Hamiltonian ([T]) we have 



In particular, for the operator function /(t = 0) = a, the evolution of /(r) corresponds to 
the evolution of a^r) = {a\a{T)\a), with the initial condition a(r = 0) = a. In this case Eq. 
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FIG. 1: Quasi-classical dynamics as described by the observable in Eq. ([5]). Parameters are e = 0.02, 
/2 = 0.01, Tn = 5, TR = 314, Tci = 2.09, |ap = 100, //ci = 1. Hence t^ < th < tr. 



dl]) can be solved exactly 



m 



a(r) = ae-^(i+^)^el"l'(^""''^-i). 



(5) 



Fig. [T] depicts the dynamics described by the observable in Eq. i^ in the coordinate- 
momentum plane. The effective coordinate is defined as x(r) = (a*(T) -|- a(r))/v^, and 
the effective momentum is defined as p(r) = i{a*{T) — a;(r))/-\/2- The corresponding clas- 
sical dynamics is described by the function aci(r) = ae~*^^"^^^'='^^, which corresponds to the 
circumference in Fig. [H Note that Eq. (j4]) maintains its form for any observable /, but 
the initial conditions for different observables are different. This is also the case for any 
quantum nonlinear Hamiltonian with many degrees of freedom. 



B. Characteristic time-scales for a closed quantum nonlinear system 

[Tl, y, y, M. In the hmit fir -^1, 



The solution ([5]) has three characteristic time-scales 
it can be re-written in the form 



a 



(r) = ad(r)e--'/2-^' [l + 0(/xr) + 0{\a\'fi'r')] . 



(6) 



The first time-scale is the characteristic classical time-scale, which can be chosen as the 
period of classical nonlinear oscillations, 

2iiuj 2ti ,^, 

^cl = = T— 1^- (7) 

The second time-scale is a characteristic time of departure of the quantum dynamics from 
the corresponding classical one 

m = ^. (8) 

This time-scale characterizes the departure of quantum dynamics from the classical one for 
classically stable systems. Historically, this time-scale was introduced for classically unstable 
(classically chaotic) systems in [15| , and it was shown to have a logarithmic dependence on 



e (see also 
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ITJI). The time-scale r?j is usually called the Ehrenfest time. The amplitudes 
of quantum and classical observables coincide at multiple times of the quantum recurrence 
time-scale, which is the third characteristic time-scale, 

TR = -. 9 

Since we are interested in the quasi-classical region of parameters, it is reasonable to 
impose the following inequalities on these three characteristic time-scales: t^ < Tf^ <^ tr. In 
our case, Td/rfi = 47r/i|a|/(l + 2fici) ~ Tr-^e ^ 1, and Tfi/r^ ^ V^/^r ^ 1- When deriving 
the first inequality, we used the conditions |ap ~ J/h = 1/e and /id — 1, which corresponds 
to the condition of strong nonlinearity. Note that the condition |a;p/i^r^ ~ 1 (see the third 
term in ([6]) in the square brackets in the expression for a(r)) gives the characteristic times 
T ^ Th, namely r/rn, = 2/e^/^ ^ 1. This means that the third term in Eq. (I6j) is small on 
the time scale r/j. For the values of parameters in Fig. [T]the inequalities r^i < t^i '^ t-^ are 
satisfied. 

C. Quantum effects as a singular perturbation to the classical solution 

As was mentioned above, the form of the differential operator K is 

K = kci + e/Xcii^q. 

The operator K^i includes only the first order derivatives and describes the classical dynamics 
of the system. Usually, the corresponding classical solution can be found by the method of 

6 



characteristics, or some alternative well-developed methods. Note that even this part of the 
solution can be rather complicated, especially for classically unstable and chaotic systems, 
and usually requires large-scale numerical simulations. (See details for closed quantum 
nonlinear systems and quantum nonlinear systems interacting with the time-periodic fields 
[7(|.) Another example which demonstrates the application of the approach based on PDEs 



with the operator K is considered in [13| for a unstable quantum nonlinear system describing 
the dynamics of a Bose- Einstein condensate with attractive interactions. 

For quantum linear systems (/id = 0) the quantum effects vanish for any values of the 
quasi-classical parameter e. The differential operator K^ includes second and higher order 
derivatives, and it describes quantum effects. The solutions of these PDEs are well behaved 
in the quasi-classical region, e <^ 1, and in contrast to the fast oscillating WKB solutions 
(typical of standard methods based on quasi-probability distributions), our method leads to 



the so-called Laplace-type expansions [10| . The crucial property of the Laplace asymptotics 
is that the dynamical observables are exponentially localized in phase space around coherent 
states. 

Quantum effects for observables represent a singular perturbation to the classical solu- 
tion. Indeed, in the quasi-classical region, quantum terms in the PDEs are represented 
by the product of the small parameter e times high order derivatives. Consequently, these 
quantum terms lead to a secular behavior of the solution, which diverges in time from the 
corresponding classical solution. Only the case e = (for finite /Xci) corresponds to the exact 
classical limit. But the problem with this limit is that for any real system e 7^ (because 
h ^ and J ^ 00). Then, even a very small value of e still "mathematically" results in a 
singular perturbation to the classical solution due to the quantum terms. 

The singularity arising from the quantum terms reminds, up to some extent, of the 
singularity provided by a "small" viscosity in the Navier-Stokes (NS) equation, describing 
the dynamics of liquid and gas flows. Indeed, in the NS equation a small viscosity multiplies 
the higher order spatial derivatives. Then, even for very large Reynolds numbers (when the 
nonlinear terms are very large compared to the viscous ones), the viscosity plays a crucial 
role in the dynamics of the flow, even though it formally represents a "small" perturbation. 
Similarly, in the quantum case the small parameter e multiplies the higher order derivatives, 
which results in a quantum singular perturbation for observables even in the "deep" quasi- 
classical region. It is this singularity that leads to a significant difference from the classical 
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FIG. 2: Frequency spectrum of the effective momentum p(t). Parameters are e = 1/900, /.t = 1/900, 
Tfi PS 15, TR « 9007r, Tci Rs 27r/3, and |qP = 900. Hence Td < r^^ ^ t-r. 

solution. 

D. Frequency Fourier spectrum for quantum observables 



The observable air) can be written in the form 



air) = Q,e~*(^+'')^"*l"'^''™(^^^) g-2|apsm2(Atr)_ 



(10) 



The first exponent in Eq. ( ITOl) is responsible for phase modulations of the classical dynamics, 
while the second one is responsible for amplitude modulations. The characteristic time-scale 
of the amplitude modulations, Tam, is defined by the condition \(y\'^fi'^Tl^ ~ 1, or by the time- 
scale Tam ~ Tfi. The time-scalc of phase modulations of the classical dynamics is defined 
by the condition |ap/irph ~ 1, or Tph ~ ^hl \f^ ^ Th- Thus, the shortest time-scale which 
characterizes the deviation of the quantum dynamics from the corresponding classical one 
is the time r^. Moreover, this time-scale is responsible for the finite width of the spectral 
line !\vn ^ 2v^/ra. 

Fig. [2] depicts the frequency Fourier spectrum of the effective momentum p(r), with 
initial condition p{0) = 0. One can see that the frequency spectrum consists of one central 
line with u = Ud = 1 + 2/ici, and a width which is approximately equal to Au ^ Auf^. In our 
case the analytical estimate gives Aun ^ 2^/2/Tt^, ^ 0.19, which is very close to the numerical 



8 



results presented in Fig. [2l Az/ ^ 0.183. The fine structure of tlie frequency spectrum is 
provided by the characteristic revival time scale tr = vr//2, or by the frequencies z/„ = 2fln, 
which are responsible for the complicated dynamics of quantum recurrences. 

III. DYNAMICS OF QUANTUM OBSERVABLES FOR OPEN QUANTUM NON- 
LINEAR SYSTEMS 

The Hamiltonian of open quantum nonlinear system interacting with an environment 
contains three terms, 

H = Hs + He + Hi^t. (H) 

The first term is typically a time-independent polynomial Hamiltonian of a general form 
which describes the self evolution of the closed system, 

-ns — 2^ J^l,sO'i • • • «Ar «i • • • ^^ , 

l,s 

where Hi^s = H^^, I = (/i, . . . , In) € Z!^ , and s = (si, . . . , sn) G Z^ . The operators ai and 
al satisfy bosonic commutation relations, [a/, a],] = 5i^k- A particular system corresponds to 
a particular choice of the coefficients Hi^s in Hs- The second term is the Hamiltonian of the 
environment, which, for example, can be modeled by a collection of harmonic oscillators, 

He = J2^<rblbg^- (12) 

q 

Usually the oscillators of the environment are assumed to be initially in thermal equilibrium, 

p,(t = 0) = Z^'e-^^/'-^, 

where Z^ = Tr[e~^^/'^^^] is the partition function of the environment, T is the temperature 
of the environment, and k^ is Boltzmann constant. The third term is the interaction Hamil- 
tonian between the system and the environment. Prototype examples are the dipole-dipole 
interaction Hamiltonian, 

Ant = ^Y1 ^".5[(«n + «n) (^ + Ml' (13) 

n,q 

and the density-density interaction Hamiltonian 

A,t = /iVa5^A,-.6tV (14) 



A. The differential operator K for many-body systems 

In a general many-body system the differential operator K can formally be written as 

Note that after explicit differentiations, exponents in K vanish. Specific examples considered 
in our previous works include: (i) a closed quantum one-dimensional nonlinear system in 
the vicinity of an elliptic [7|, [ll| or a hyperbolic ll|, [l3| point; (ii) chaotic systems describing 



the interaction of atoms with radiation and external radio frequency fields [7[ ; and (iii) the 



ields [71; 



quantum Brownian motion problem for a nonlinear system oscillator 



a. 



B. Frequency Fourier spectrum of p(r) in the presence of an environment 

Let us introduce formally a relaxation (dissipation) term into Eq. ^. Namely, we con- 
sider the function 

«(r) = a e-^--*(i+'^)- ^|„P(e-2..._i)^ ^^^^ 

where the parameter 7 plays the role of an effective relaxation. The characteristic time scale 
of relaxation is r^ = I/7. We consider the frequency Fourier spectrum of the momentum 
p{t) = i{a*{T) — a(r))/-\/2, with a(r) given by Eq. ( fT6l) . for two cases: (i) Tj ^ th (Fig. [3k), 
and (ii) r-^ < r^ (Fig. [3b) (similar dependencies can be built for the effective coordinate 
x(r)). As one can see, when the influence of the effective dissipation is small (Fig. [3^), the 
width of the Gaussian spectral line (at the level e~^) is still determined by the time-scale Tn 
{Aufi, ~ 2-\/2/Tri ~ 0.19), and not by the environment (Ai^^ ~ 27 = 0.001). The numerical 
results give Au ^ 0.186. Note that in this case the fine structure of the spectral line is not 
completely destroyed, as both time-scales, tr ~ 2826 and r^ = 2000, are of the same order. 
In the case of strong dissipation (Fig. [3b), the width of the spectral line has a Lorentzian 
form, 

Re(p,)=72Re(po)/(7' + i^'), 

with a width (at Re(pj,) = 1/2) determined by the dissipation parameter 7 (Az/^ :^ 27 = 1). 
The numerical results are in good agreement, Au ^ 1. Also, the fine structure is destroyed, 

10 



as in this case r^ = 2 ^ tr ~ 2826. In |8| we studied the concrete example of the QNO 
interacting with an environment in which both relaxation and decoherence take place, and 
we found that the frequency spectrum behaves in a similar way as the toy model discussed 
in this subsection. 

C. Characteristic parameters for observation of quantum effects after decoherence 
and relaxation 

As was discussed above, for the simple closed quantum nonlinear system given by Eq. ([1]) 



there are three characteristic time-scales (see 



10| for details on multi-dimensional systems). 



Due to the interaction with the environment, two new time-scales appear: Td - a very short 
decoherence time, and r^ -the relaxation time. All of these five time-scales depend on the 
parameters of the system and the environment. The typical region of parameters in which 
one can observe quantum effects after decoherence and relaxation is ra ^ Td < r^ < r^ < 
tr. In the following we will consider a system which satisfies this region of parameter 
a "quasi-classical nonlinear mesoscopic system". The key inequality is Tfi < r^. In this 
case, the deviation of the quantum dynamics from the classical one formally works as an 
effective "quantum relaxation" (or a "quantum amplitude modulation"), which gives the 
main contribution to the frequency spectral line width. The relations between Tc\ and Tf^ , 
and between tr and r^ are not so important. There can be additional time-scales related 



to accumulation of quantum phases 



multi-dimensionality 10|, etc. The details for a 



one-dimensional case were presented in |8|, |9|| . 



IV. AN EXACT SOLVABLE EXAMPLE OF AN OPEN QUANTUM NONLINEAR 

SYSTEM 

A. Phase decoherence 

Although the PDEs described above look rather complicated, especially for open quantum 
nonlinear systems, we have found the exact solution for a quantum nonlinear oscillator 
interacting with the environment in the special case of a density-density type of interaction, 
as in Eq. (ll4l) . The type of interaction does not provide relaxation processes through energy 
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FIG. 3: Fourier frequency spectrum of the momentum p{t) obtained from Eq. ()16p . Parameters 
are: a) 7 = 0.0005, r^ = 2000 ^ r;^ « 15; b) 7 = 0.5, r^ = 2 < r/j ss 15; all other parameters are 
the same as in Fig. [2j 

exchange between the QNO and the environment, but leads to phase decoherence. These 
effects result in the decay of the amplitude of oscillations of the QNO (similar to the effects 
of relaxation), and survive in the classical limit, where they correspond to the dephasing of 
the QNO. 

We summarize here the results of 9j in the context of the quantum-classical transition for 
observables and the frequency Fourier spectrum. We choose the system Hamiltonian Hs as 
in Eq. ([1]), the environment Hamiltonian H^ as in Eq. (fT2l) . and a density-density interaction 
Hamiltonian i^^int as in Eq. flMl) . For the model under consideration, the interaction with 
the environment introduces a single time scale, ra , which plays the role of a decoherence 



time. This time-scale is not small, and survives even in the classical limit: |a| 



00, 



\f3. 



00, ^ — > 0, ^|ap = J = const, hlP^l"^ = J^= const, X^ = const. The typical region 



of parameters in which one can observe quantum effects is Tc\ < r^i < t^ < tr. The key 
inequality is now r^^ < r^. In this case, the deviation of the quantum dynamics from the 
classical one formally works as an effective "quantum relaxation" (or a "quantum amplitude 
modulation"), which gives the main contribution to the frequency spectral line width. The 
relations between Td and t^ , and between tr and Td are not so important. 

Although this model is rather trivial because the full Hamiltonian can be diagonalized in 
the number basis for the joint system-environment Hilbert space, it is a useful model system 



12 



for the purposes of demonstration of our approach. Following our previous results 91] it is 
possible to write an exact linear PDE for any quantum djTiamical observable in the joint 
Hilbert space 

f{a\a-pt,(3,^-t) = {a,P,^\f{t)\a,(3^), (17) 

where f{t) = f{a^{t),a{t);bjt),bq-{t)) is a generic Heisenberg operator function, and 
\a, jSq) = \a) Ylj IPq^) is an initial coherent state of the system and the environment. Here 
a^(t), a{t), bJt), and bg-(t) are the Heisenberg bosonic creation and annihilation operators 
for the system and the environment, respectively, and q = {qi,q2, ■ ■ •)• The corresponding 
PDE has the form 



-/(a*,a;/5|,/3,-.) =i^/(a*,a;/5|,/5,-;t), (18) 



d_ 

Wr 

where the differential operator K includes the derivatives of different orders over a*, a, 
/3| and (3^, and depends on the explicit form of the corresponding full Hamiltonian. As 
before, the general form of the differential operator K i^ K = K^x + iTq. The operator Kc\ 
includes only the first order derivatives and describes the classical dynamics of the system 
and environment. The operator K^ describes the quantum effects of the system and the 
environment. The explicit expressions for both these operators are given in 9|. 

In order to study the reduced dynamics of the system, the function /(a*, a; /5i, /Jg-) has 
to be traced over the variables of the environment /3i, (3^. We have assumed above that 
initially each environmental oscillator is populated initially in the coherent state |/3q>). Let 
us now assume that the each environmental oscillator is initially in a mixed thermal state 
at temperature T. Then we should perform an additional averaging of the environmen- 
tal oscillators over the thermal distribution. The corresponding procedure is thoroughly 
explained in |9|. The exact solution for the system observable (Q;(r))f, averaged over the 
environmental variables, is 

{a{r))s = air)RiT), (19) 

where a(r) is defined in Eq. (3), and 

R{t) = 1[r'-'^{t), 

Q 

R'-'^i'r) = , ,, ^ .,, ; ■ (20) 
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B. Classical limit 

Let us write the complex quantity R^^ (r) in terms of its modulus and phase, R^^ (r) = 
e^'/'<'^M|/2('?1(r)|. Then, we have from Eq.(|2Ul) 

\R{r)\ = e-^M, 

V r°^ 

r(r) = -—j^ dqq'H\R^'^\r)\), (21) 

where V is the volume of the thermal bath. In the classical limit {h -^ 0) we have for r(r) 



and the phase if^^ij) is 



2rr rj-2n^ u^ J, -u;. 



Jq)( -\ - , ^qf^BT 



dq^, (22) 



<^(«)(^)~_:x:i^r. (23) 

UJq'UJ 

The function a(r) in Eq. (fT9l) coincides with that in Eq. ([5]). It is clear from Eq. (TT9|) 
that under the condition 

r/i < Trf, (24) 

the width of the frequency spectrum of (^(r))^: is defined by the time-scale, r^, and not by 
the interaction with the environment. In the opposite case, r?j ^ r^, the width of the spectral 
ine is determined by the interaction with the environment. A similar result was obtained in 
for the QNO interacting via the dipole-dipole interaction with the environment Eq. (TT3!l . 
But in the latter case, the time-scale r^ in Eq. (I2ll) should be substituted by the relaxation 
time r-y. 

V. ESTIMATES FOR CONCRETE SYSTEMS 

Our main statement is that generally there is no classical limit for the dynamics of 
quantum nonlinear systems interacting with the environment, even when these systems are 
in the deep quasi-classical region of parameters. The corresponding systems were called 
above quasi-classical nonlinear mesoscopic systems (QCNMS). In this context we note that 
most classical systems surrounding us represent a very particular exception due to (i) either 
an extremely deep quasi-classicality (extremely small value of e) and/or (ii) a very strong 
interaction with the environment. At the same time, the general belief in the recent scientific 
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literature is that after the process of decoherence, the quasi-classical system can be described 
by using classical probabilistic approaches. According to the results discussed here it appears 
to be true only (i) for quantum linear systems (with quadratic Hamiltonians) or (ii) for 
quantum nonlinear systems with significantly small value of a quasi-classical parameter e. 
For the QCNMS quantum effects survive after the processes of decoherence and relaxation 
took place. Moreover, these quantum effects make a crucial contribution to the dynamics 
of observables. This observation may have significant relevance for the understanding of 
the properties of noise in complex quantum systems and nanodevices. In particular, the 
performance of future BEC based interferometers and nano machines will be limited by the 
level of noise. 

The key condition for survival of quantum effects for observables related to the time-scale 
Th is Th < T^i which, in the simplest case of the quantum nonlinear oscillator can be written 
in the form 

e = ^ = 2/ieiv^r^ > 1. (25) 

We now present estimates for different real QCNMS that may satisfy the above condition, 
and therefore may lead to the observation of certain quantum effects that survive the process 
of environment- induced decoherence and dissipation. 

A. Bose-Einstein condensates in a one-dimensional toroidal geometry 

We start with a one-dimensional BEC confined in a toroidal geometry, and described by 
the quantum field equation (see (l3|, [l8|, and references therein) 



dr 



o2 

90-2 + 2^^^^^ 



^. (26) 



Here e = 4Ra/S, R is the radius of the toroidal trap, S is the area of the cross-section of the 
torus, and a the interatomic s-wave scattering length (a > for a repulsive interaction, and 
a < for an attractive interaction). The dimensionless time is r = ht/2mR^ . The operator 
\l/(6', r) can be expanded as 



^ie.T) 



^ CO 

= Y. a,(r)e'^^ (27) 



2^i. 

fe=— oo 
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Here ctkir) and alir) are annihilation and creation bosonic operators, respectively, and the 
field operator is periodic \l/(^ + 27r, r) = \l/(^, r) and satisfies the normalization condition 



/ ^\9,T)^i9,T)d9= J2 ^k = N, (28) 

■^0 fc=-oo 

where fik is the operator of the number of particles in the mode with momentum k, and A*" 
is the operator of the total number of particles. 

In the following we only consider the case of repulsive interactions, a > 0. /,From Eqs. fl2^ 
and fl27j) it follows that the operators dk{T) satisfy the following system of coupled first-order 
differential equations: 

oo 

idk = k dk + e 2_^ (l'ki^k2(^k-sSk+ki-k2-k3,0, (29) 

ki,k2,k3=—oo 

where "dot" means derivative with respect to r. Here we shall limit ourselves to consider 



only a single mode in Eq.( 
a more general case see (|l9l. 



d'k = dk5k',ki which is stable under the condition a > (For 



20| ) ) ■ In this simplified case Eq. fl29l) takes the form 
idk = k'^dk + edldl = [dk, Hcs], (30) 



.2 ^ UtA I ^ {?A?, \2 



with the effective Hamiltonian 

^cff = ik'^ - ^jdldk + ^{dldkY- (31) 

To solve the system Eqs. (l30|) . (!3T]) we use the above described techniques of projection 
onto the basis of coherent states. Let us assume that at r = the kth mode of the bosonic 
field can be represented by a coherent state, \ak), described by a complex number a^. We 
denote 

ak{r) = {ak\ afc(r) \ak) = ak(t, ak, a^*). (32) 

Note that all atoms A^ occupy the single mode k, that is {ak\nk\<yk) = nk = N. The exact 
linear PDE for the observable afc(r) is 

ak{r)=Kak{T), 
akiO) = ak, 



where 

d 



K = i ( /c^ + 5|afcp ) ( a*k-K^ - cc. 



■)(■ 



i2 



■2V^"- d{aiy 



+4f(«^)'^---)- (34) 
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It is more convenient to write these equations using action-angle variables. Namely, instead 
of the variables Uk and a^ we use the variables Uk (remember that in this simplified case the 
number of atoms in mode k is fixed, n^ = A^) and Ok-, where 



«fc = v^e-'"^ (35) 

Using the expressions 

( * 9 \ . d 

92 \ 9 ., 92 



'at? ^/ ., - c.c. ] = i-— - 2iN 



one can derive the following equation for ^^(r) in new variables 

Or V 2 y d9k ON 89 ^ ' 

This equation possesses a solution of the form of a finite amplitude periodic wave 

akir) = exp | - ik'^r - (1 - exp{-ieT))\ak\'^ jak- (37) 

This solution has two characteristic time-scales 

1 27r 

^1 = 1 — fTT' '^R = n- (38) 

a/c k \e\ 



The first one describes the breakdown of quantum-classical correspondence, and the second 
one is the time-scale of quantum revivals |12l . Il3l . Il4l ]. Note that Eq. fl37|) formally turns 



into the GP solution (which we also will call a "classical" field theory solution) 

a^'(r) = exp I - i(k'^ + e\ak\nT\ak, (39) 

when \e\ — > 0, |afcp = A^ ^ 00, and |e||afcp = const. 

For this one-dimensional BEG system the condition Eq. (!25|) for observation of quantum 
effects after decoherence and relaxation is reduced to the following 

mRS , , . 

tn = T^TTT- « ^7. 40 

2n\a\a 

where t^ is the relaxation time, and m is the mass of the BEG atom. To estimate the time- 
scale tfi we assume that A^ = 10^ ^^Rb atoms (a = 2.5 x 10~^cm) are trapped in a toroidal 
trap with radius R = 5 x 10~^cm and cross-section S = 10"'^cm2, which implies tfi ~ 4.5ms. 
The corresponding bandwidth of the frequency spectrum, which characterizes the quantum 
effects related to the time-scale tn, is Au ^ 2\/2/tn, ~ 0.6kHz. 
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B. Relation between the time-scale t/j and phase diffusion of two Bose-Einstein 
condensates 

The effect of phase diffusion of the relative phase between two BECs due to atomic 



coUisions was studied theoretically in |21| and was rece ntly observed in the high atomic 



density regime with two BECs trapped on an atom chip 221 ] . 

Let us summarize here the main ideas behind phase diffusion in the simplest ideal case. 
Imagine a Bose-Einstein condensate that is symmetrically split into two pieces via a double- 
well potential. Assuming that the split process is slow enough (i.e., the barrier is raised on a 
time scale long compared to the inverses of the excitation frequencies of the initial potential 
well), but fast enough to freeze the relative phase between the two condensates in each well, 
the final state of the condensates after the split can be described as a state \ip) that is a 
superposition over many relative number states 

^ / AM 

where N is the total number of atoms and ip is the relative phase between the condensates 
in each well. Here we have assumed that during the split process the atomic interactions 
are negligible. For simplicity we will also assume that the relative phase between the two 
condensates is zero, c/) = 0. After the split, each condensate evolves independently (the 
barrier is sufficiently raised to suppress tunneling between the wells). Because of atom-atom 
interactions, the energy of number states E{k, N — k) have a quadratic dependence on the 
atom numbers in each well, k and N — k, so that the different relative number states have 
different phase evolution rates. The state vector (HTi) evolves as 



where u = 2E{N/2)/h is the frequency of each well, and ^ = \ ^^2 |fc=JV/2 = 2/i^ is the 

effect of nonlinearities. To study the phase distribution of the evolved state one projects 

this evolved state onto phase states. These are orthonormal states of the form 

1 ^ 

\<P,) = ^==Y^e^f^%N-k), (43) 

with (j)p = 2-Kp/{N+l) and p = -N/2, . . . , A^/2. In the limit iV > 1, the phase distribution of 



the state (gj]) is P(0) = |(0|x,t)|2 = ^/7rJ2{A(j)f ex.p{~(j)^/2{A(f)f), with a phase dispersion 
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that evolves in time as 

Acj>{tf = A<Pl + RH\ (44) 

where A</)q = 1/A^ is the phase dispersion for the initial two-model coherent state pTj) . and 
R = vN^ is the rate of phase diffusion. This rate defines a phase diffusion time-scale 

^ph.diff = r— , (45) 

2/ivA' 
which coincides with the time-scale r^^ in Eq.Q of breakdown of quantum-classical corre- 
spondence of the quantum nonlinear oscillator initially prepared in a coherent state with 
mean number of excitations A^ = |ap. 

We conclude from the above considerations that an alternative way to observe the effect 
of the "quantum" time scale r^ in the dynamics of the quantum nonlinear oscillator (QNO) 
is to analyze phase diffusion of two condensates (which can be modeled as two uncoupled 
QNOs after the splitting process), initially prepared in a quasi-classical coherent state, as 



in Eq. 04ip . A related experiment was performed in [22| for an initial two-mode number- 
squeezed state instead of a two-mode coherent state. In that case, the initial phase dispersion 
is much wider, A</)q ~ s/N, and the rate of phase diffusion is much larger, R ^ sy/N^, where 
s ~ viV ^ 1 is the squeezing parameter. 

C. The time-scale Tfi for mechanical resonators and for nonlinear optical systems 

For a mechanical resonator or cantilever the quasi-classical parameter is e = 1/n, where 
n is the average number of levels involved in the quantum state of the resonator, that we 
assume to be a coherent state. The dimensionless relaxation time is T-y = 2Q, where Q is 
the resonator's quality factor. The condition Eq. fl25|) takes the form 



Ocantilever = ^ » 1- (46) 



Different aspects of cantilevers, from kilohertz to gi gah ertz frequencies, including their non 



24| 



linear properties, are discussed, for example, in [23|, 

A condition similar to Eq. (1461) holds for quantum nonlinear optical systems in high 
quality resonators. In this case, n is the average number of photons in the initially coherent 
state of the cavity resonance mode, and the classical parameter of nonlinearity can be written 
as fid = xJjj^ciw, where x is the nonlinear susceptibility, and cUcav is the cavity resonance 



frequency [25 1. 
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VI. CONCLUSIONS 

In this paper we have reviewed the effects of singular perturbations resulting from quan- 
tum terms in the dynamical equations for observables of open quantum nonlinear quasi- 
classical systems. We have argued that when the time-scale for quantum-classical depar- 
ture, generally given by the time-scale r/j, is much shorter than the dissipation time scale 
Xy, certain quantum effects survive the process of decoherence, and could be observed from 
characteristic properties of the time-evolution of observables, such as in the frequency spec- 
trum and in the noise spectrum. With recent advances in quantum technology we expect 
that the key condition (1251) for detecting such effects may be experimentally realized, and 
quantum effects related to the time-scale Tn, can be observed in the quasi-classical region of 
parameters. 
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